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BIOLOGICAL DATA 

This invention relates to the compression of biological sequence data for electronic storage. 

The nature of biological sequence data (eg. DNA and protein sequences) means that electronic 
storage is perfectly suited. Not only does the sheer volume of data necessitate large-scale 
5 storage, but electronic storage allows rapid and efficient searching of the data eg. for 
homologous sequences. Since the advent of initiatives such as whole genome sequencing, the 
amount of storage required has increased significantly. The storage requirement for the yeast 
genome, for instance, is huge. Whilst large capacity storage systems continue to fall in price, 
one of the rate-limiting steps when dealing with sequence data is the transfer from storage 
10 medium into memory (eg. hard drive into RAM) and any developments which significantly 
reduce the size of sequence data files would be welcomed. 

Biological sequence data is typically . represented in an alphabetic manner, rather than by 
chemical formula, with each letter representing a monomer unit in a biological polymer (Table I). 
For instance, DNA sequences are represented as strings of letters chosen from a simple 
15 four-letter "alphabet". Each A, C, G or T represents a monomer unit (nucleotide) in a DNA 
polymer. Similarly, proteins are made up of twenty different monomer units (amino acids), 
which have each been assigned single letter codes. 

Because of its alphabetic nature, biological sequence data is naturally suited to electronic storage 
in alphabetic text form. Alphabetic text-based computer information is generally stored and 
20 manipulated using the char datatype, using 8 bits (1 byte) and a conventional file of biological 
sequence data is made up of a string of characters of datatype char. A conventional file of 
sequence data uses a single byte to represent each monomer, so the amino acid sequence of the 
glycogen synthase protein, for example, requires 737 bytes of storage using the one-letter amino 
acid code, and the corresponding DNA sequence requires 221 1 bytes. 

25 The char datatype, however, was designed for representing a full character set, including upper 
and lower case letters plus numbers, punctuation, and other characters, and each 8-bit char can 
represent 256 (2 8 ) different values. Using the char datatype and alphanumeric characters to store 
DNA sequences therefore fails to utilise 252 of the available values. Similarly, protein sequences 
waste 236 values. Other datatypes which are in common usage for data storage include int (16 

30 bits), long (32 bits), float (32 bits), although this may vary from machine to machine. 



<WO 9901940A1_I_> 



WO 99/01940 



-2- 



PCT/GB98/01937 



The DNA and RNA alphabets each consist of 4 letters and, rather than storing these sequences in 
alphanumeric form using strings of char datatype, using a sub-byte datatype would enable a 
significant storage saving. Degenerate nucleic acid sequence information (which can be 
represented using a 16 letter alphabet) and protein sequences could also be treated in this way. It 
5 would therefore be useful to define a sub-byte datatype in order to take advantage of the small 
size of the biological alphabet. 

The commonly used blast sequence comparison program converts single byte char data into a 
half-byte working space whilst manipulating data. This is a temporary measure, however, and 
data is not stored in this manner using a specific sub-byte datatype. 

10 The invention is based upon the realisation that using a whole byte to represent a monomer in a 
biological sequence is not the most efficient means of permanent storage. 

According to the invention, there is defined a sub-byte datatype for the storage or manipulation 
of biological sequence data in a programming language or a database. 

The invention also provides a programming language or a database which utilises a sub-byte 
15 datatype for the storage or manipulation of biological sequence data. 

According to a further aspect of the invention, there is provided the use of a sub-byte datatype in 
the storage or manipulation of biological sequence data. 

By "sub-byte" it is meant fewer than 8 bits. 

The datatype may be intrinsic to a program or programming language, or it may be user-defined. 
20 The invention is not limited, however, to situations where a formal datatype must be defined. 

According to a further aspect of the invention, there is provided a computer program which 
stores biological sequence data using fewer than 8 bits to represent each monomer in said 
sequence data. 

The invention also provides a file containing biological sequence data, wherein each monomer in 
25 said sequence data is represented using fewer than 8 bits. 

According to a further aspect of the invention, there is provided a method for compressing 
biological sequence data, comprising representing each monomer in said sequence data by using 
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fewer than 8 bits. 

The invention also provides a method for reducing the size of a file in which biological sequence 
data is represented using 8 or more bits per monomer, comprising replacing the representation of 
each monomer with a representation using fewer than 8 bits. 

5 According to a further aspect of the invention, there is provided a computer programmed to store 
biological sequence data by using fewer than 8 bits to represent each monomer in said sequence 
data. 

According to a further aspect of the invention, there is provided a computer comprising means 
for alphabetic entry of biological sequence data, means to convert said sequence data into a 
10 format wherein each monomer unit is represented using fewer than 8 bits and, preferably, means 
to store said data. 

According to a further aspect of the invention, there is provided a storage medium holding 
biological sequence data, wherein said sequence data is stored using fewer than 8 bits to 
represent each monomer in said sequence data. 

15 The storage medium may be in any appropriate form, such as a floppy disk, a CD-ROM, or a 
fixed disk drive. 

According to a further aspect of the invention, there is provided a method for transmitting 
biological sequence data, comprising compressing the data by representing each monomer in said 
sequence data by using fewer than 8 bits before transmission, for instance over a network. 

20 According to a further aspect of the invention, there is provided biological sequence data which 
has been electronically stored using less than 8 bits to represent each monomer in said sequence 
data. 

The biological sequence data may be of any suitable kind, such as DNA sequence, RNA 
sequence, and protein or polypeptide sequence. 

25 It will be apparent that nucleic acid sequences can be represented using 2 bits to represent each 
monomer (nucleotides A, C, G, or T/U). Accordingly, a 2 bit datatype may be defined according 
to the invention for the storage or manipulation of nucleic acid sequences. Such a datatype is 
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referred to herein as base. 

By representing each nucleotide in a nucleic acid sequence by using only 2 bits, 4 nucleotides 
can be stored in a single byte. This represents a 75% compression compared with the 
conventional representation of each nucleotide using a single byte. 

5 Where a nucleic acid sequence is not definite, more than 2 bits are required to represent each 
nucleotide. For instance, where a nucleotide has not been unequivocally determined, the symbol 
"N" is used according to IUPAC convention. The alphabet of this IUPAC convention (Table I) 
has 16 members. This can be conveniently represented using 4 bits per member. Accordingly, a 4 
bit datatype may be defined according to the invention for the storage or manipulation of 
10 degenerate or uncertain nucleic acid sequences. Such a datatype is referred to herein as longbase. 

By representing each nucleotide in a sequence by using 4 bits, 2 nucleotides can be stored in a 
single byte. This represents a 50% compression compared with the conventional representation 
of each nucleotide using a single byte. 

As an alternative to using 4 bits to represent degenerate or uncertain nucleic acid sequences, 
15 under certain circumstances these features may be accommodated where 2 bits are used, as in 
base. For instance, where a DNA sequence is stored in a data file using 2 bits per nucleotide, 
parallel files could be utilised which contain "modifying" data to qualify details in the sequence 
file. For instance, the second file may contain an indication that whilst nucleotide 221 is given as 
guanine in the sequence file, in fact it may be any purine. Obviously, the choice of using such a 
20 "modifying" file or using more than 2 bits to represent the sequence depends on the particular 
situation, but the choice is routine. 

It will further be apparent that protein sequences require at least 5 bits to represent each monomer 
(20 amino acids) since 2 4 =16 and 2 5 =32. Whilst this is encompassed within the invention, 5 bits 
is an awkward length, being an odd number. 6 bits is more convenient and, furthermore, this 
25 allows a degree of degeneracy to be incorporated into the sequence (2 6 =64). Accordingly, a 6 bit 
datatype may be defined according to the invention for the storage or manipulation of protein 
sequences. Such a datatype is referred to herein as aminoacid. 

By representing each amino acid in a protein sequence by using 6 bits, 4 amino acids can be 
stored in 3 bytes. This represents a 25% compression compared with the conventional 
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representation of each amino acid using a single byte. 

The degree of degeneracy incorporated into a 6-bit representation or datatype also allows an 
amino acid to be represented in terms of codons, of which there are 64. A datatype used in this 
way is referred to herein as codon. Each single codon value represents a single codon, which 
5 inherently also defines an amino acid. In effect, the codon datatype represents three base entries, 
just as a codon is made up of three nucleotides. By using 6 bits to represent each codon, 4 codons 
can be represented in 3 bytes. This represents a 75% compression compared with the 
conventional representation of each codon using 3 bytes. It will also be appreciated that a full 
byte could be used to represent each codon, which would allow a degree of degeneracy and 
10 would represent a 67% compression compared with using 3 bytes to represent each codon. 

It should be bome in mind that the various datatypes and compressions described above may not 
be suitable in all circumstances. For example, the programming language C requires a string to 
have a NULL terminator. This is not possible with the base datatype, for instance, because all of 
the 4 possible values (permutations of 2 bits) are used to represent information, which does not 
15 allow a terminator to be represented. 

Similar caveats apply to longbase. The IUPAC convention uses 15 representations for a DNA or 
RNA sequence, which does allow the sixteenth permutation to represent a terminator. In certain 
circumstances, however, a value may be needed to represent a gap (representing an unknown 
sequence of unknown length) which would remove the possibility of having a terminator. The 
20 codon datatype is also "full" since each of the 64 available values represents a codon. 

Whilst these datatypes may not be universally applicable, however, they are not without utility 
since not all programming languages or databases have such a terminator requirement. A further 
problem in using the datatypes of the invention in languages such as C is the international ANSI 
standard which does not recognise these datatypes. However, new languages, such as Java which 
25 is still in early development, currently have less strict standards and may be amenable to the 
introduction of new datatypes at this stage. 
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TABLEI 

IUB/IUPAC standard biological sequence codes 

Single letter nucleotide codes 

A Adenine C Cytosine 

5 G Guanine T Thymine 

U Uracil 
Degenerate nucleotide codes 

In addition to the five above codes: 



N any (A/C/G/T) 


R 


puRine (G/A) 


Y 


pYrimidine (T/C) 


K Keto (G/T) 


M 


aMino (A/C) 


S 


Strong (G/C) 


W Weak (ATT) 


B 


not A (C/G/T) 


D 


not C (A/G/T) 


H not G (A/C/T) 


V 


notT(A/C/G) 






Single letter amino acid 


codes 








A Alanine 


C 


Cysteine 


D 


Aspartate 


E Glutamate 


F 


Phenylalanine 


G 


Glycine 


H Histidine 


I 


Isoleucine 


K 


Lysine 


L Leucine 


M 


Methionine 


N 


Asparagine 


p Proline 


Q 


Glutamine 


R 


Arginine 


S Serine 


T 


Threonine 


V 


Valine 


W Tryptophan 


Y 


Tyrosine 







In addition: 

B represents asparagine or aspartate ie. N or D 
Z represents glutamine or glutamate ie. Q or E 
U represents selenocysteine 
25 X represents "any amino acid" or "unknown" 
* represents a translation stop 
- represents a gap of indeterminate length 
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CLAIMS 

1. A sub-byte datatype for the storage or manipulation of biological sequence data in a 
programming language or a database. 

2. A programming language or a database which utilises a sub-byte datatype for the storage or 
5 manipulation of biological sequence data. 

3. The use of a sub-byte datatype in the storage or manipulation of biological sequence data. 

4. A file containing biological sequence data, wherein each monomer in said sequence data is 
represented using fewer than 8 bits. 

5. A method for compressing biological sequence data, comprising representing each monomer 
10 in said sequence data by using fewer than 8 bits. 

6. A method for reducing the size of a file in which biological sequence data is represented using 
8 or more bits per monomer, comprising replacing the representation of each monomer with a 
representation using fewer than 8 bits. 

7. A computer programmed to store biological sequence data by using fewer than 8 bits to 
15 represent each monomer in said sequence data. 

8. A storage medium holding biological sequence data, wherein said sequence data is stored 
using fewer than 8 bits to represent each monomer in said sequence data. 

9. Biological sequence data which has been electronically stored using less than 8 bits to 
represent each monomer in said sequence data. 
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